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Nuclear quadrupole moments (Qs) in three isotopes of potassium (K) with atomic mass numbers 
39, 40 and 41 are evaluated more precisely in this work. The Q value of 39 K is determined to 
be 0.0614(6) 6 by combining the available experimental result of the electric quadrupole hyperfine 
structure constant (B) with our calculated B/Q result of its 4P 3 / 2 state. Furthermore combining 
this Q value with the measured ratios Q( 40 K)/Q( 39 K) and Q( 41 K)/Q( 39 K), we obtain Q( 40 K) = 
-0.0764(10) b and Q( 41 K) = 0.0747(10) b, respectively. These results disagree with the recently 
quoted standard values in the nuclear data table within the given uncertainties. The calculations 
are carried out by employing the relativistic coupled-cluster theory at the singles, doubles and 
involving important valence triples approximation. The accuracies of the calculated B/Q results 
can be viewed on the basis of comparison between our calculated magnetic dipole hyperfine structure 
constants {As) with their corresponding measurements for many low-lying states. Both A and B 
results in few more excited states are presented for the first time. 

PACS numbers: 21.10.Ky, 31.15.aj, 31.30.Gs, 32.10.Fn 



I. INTRODUCTION 

Potassium (K) atom has three naturally abundant iso- 
topes with atomic mass numbers 39, 40 and 41. Using 
the modern femtosecond laser frequency combs and po- 
larization quantum-beat techniques, high precision mea- 
surements of hyperfine structure constants in 4P and 3D 
states are carried out Hi- Also, a number of measure- 
ments of these quantities were carried out in the ground 
and other states long ago using the atomic beam mag- 
netic resonance and level crossingtechniques (e.g. see 
review article by Arimondo et ai. Q). Theoretical stud- 
ies of these quantities are of immense interest to atomic 
physicists to test the accuracies of the wave functions 
in the nuclear region [4j-|6|. However, theoretical evalu- 
ation of these quantities require atomic calculations and 
nuclear moments 0, 0, 0, Q • Nuclear magnetic moments 
(/is) of the above K isotopes are known very precisely and 
the reported results from various studies matches reason- 
ably well with each other Q. On the other hand, the 
reported nuclear quadrupole moments (Qs) from various 
works on these isotopes differ significantly. For example, 
Q value of 39 K is reported as 0.07(2) b [n|, 0.049(4) b 0, 
0.0601(15) b [ll| and 0.0585 b within one percent error 
The latest result, 0.0585(6) b, is now considered as 
the standard Q value for 39 K Accurate knowledge 

of Q values of these isotopes are useful in many applica- 
tions. These information are interesting in order to test 
the potential of nuclear models 14|, HH , acquiring infor- 
mation about local symmetry [16| , to hnd out asymmetry 
parameters in nuclei [13, EH , for studying the Mossbauer 
spectroscopy for the structural determination of the ele- 
ment containing solid state compounds [l9[ etc. 

In this paper, we analyze the electric quadrupole hy- 



perfine structure constants for many states in K and re- 
port precise Q values of its above mentioned isotopes. As 
discussed later, we find the new values to be larger than 
the considered standard values in the literature. Atomic 
wave functions are calculated using relativistic coupled- 
cluster (RCC) method in the Fock-space representation 
and matrix elements of the hyperfine interaction Hamil- 
tonians are estimated using these wave functions in the 
considered atom. 

The rest of the paper is organized as follows: In the 
next section, we present briefly the theory of hyperfine 
structure in an atomic system and the single particle ma- 
trix elements of the interaction Hamiltonians which are 
used to evaluate the hyperfine structure constants. In Sec 
III, we explain the RCC method little elaborately for the 
calculation of atomic wave functions. Then, we present 
the results and their discussions before summarizing the 
work. Unless stated otherwise, we use atomic unit (au) 
throughout this paper. 



II. THEORY OF HYPERFINE STRUCTURE 

The hyperfine structures of energy levels in an atom 
arise due to the interaction between electron angular mo- 
menta with the nuclear spin. Details of this theory is 
given by C. Schwartz in a classic paper [20(. Mathemati- 
cally, the hyperfine interaction Hamiltonian is given in a 
general form as non-central interaction between electrons 
and the nucleus in terms of tensor operators as 



(2.1) 
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where T e and T„ arc the spherical tensor operators 
of rank k in the space of electronic and nuclear coordi- 
nates, respectively. In the first order perturbation theory, 
the hyperfine interaction energy Wf of hyperfine state 
\F; I J) with total angular momentum F = I + J for / 
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and J being the nuclear spin and electronic angular mo- 
mentum of the associated fine structure state \J,Mj), 
respectively, taking up to k — 2 is given by 



W f = ^AR 



\R(R + 1) - 21(1 + 1) J( J + 1) 



B2 21(21 -1)2 J (2 J -1) ' (2 ' 2) 

with R = F(F + 1) - 1(1 + 1) - J( J + 1), and A and B 
are known as the magnetic dipole and electric quadrupole 
hyperfine structure constant for k = 1 and k = 2, respec- 
tively. The advantage of expressing the change in energy 
in this form is it separates out the electronic and nuclear 
factors for which the calculations can be carried out in a 
simple approach. Here A and B are given by [2(J HH 



A = hn9i 



(J\\TP\\J) 



and 

B = 



Q 



V / J(J+1)(2J+1) 
8J(2J-1) 



(2.3) 



(2J + l)(2J + 2)(2J + 3) 



(J\\TP\\J). (2.4) 



In the above expressions, /in and gi = fi/I are the 
nuclear magneton and gyromagnetic ratio, respectively. 
Since our intention is to verify accuracies of Q values, we 
estimate B/Q results in this work. 

The reduced matrix elements of the electronic spher- 



ical tensor operators, = ^tt 1 , in terms of single 
orbitals are given by [2(], Hl[ 

(«/ll«i 13 l^} = -(«/ + ^)(-«/l|c (1) ll^> 

I 



(fc) 



~*e&+9* (2.5) 



and 



dr 



PfQi + Q/Pi 



where Ki and Pi (Qi) are the relativistic angular momen- 
tum quantum number and large (small) component of 
Dirac spinor for the corresponding orbital i, respectively. 
The reduced matrix elements of Racah tensors (C^) are 
given by 22] 



(KfWC^Wm) = (-i)^ +1 / 2 v /(2 J/ + i)(2 j4 + i) 

jf k ji 



1/2 1/2 H'/'*' ^ 



with the angular momentum selection rule n(£f,k,£i) 
1 when if + k + £i = even for the orbital angular momen- 
tum if and li\ otherwise it is zero. 



III. METHODS FOR CALCULATIONS 

A. Single particle orbital generation 

Accurate generation of atomic orbitals in the nuclear 
region is very important for the present study. We 



TABLE I: Used £ and v parameters for different symme- 
tries to construct GTOs. 





s 


P 


d 


f 


g 


c 


0.0002 


0.0004 


0.0003 


0.0005 


0.0004 


V 


1.917 


1.79 


1.77 


1.76 


1.75 



consider here Gaussian type of orbitals (GTOs) which 
provide natural description of relativistic wave-functions 
within nucleus [23|-[25| as basis to construct the mean- 
field orbitals in the Dirac(Hartree)-Fock (DF) approach. 
Kinetic balanced condition between the large and small 
components of Dirac spinor are imposed to ensure correct 
non- relativistic behavior of the orbitals [25, 26]. GTOs to 
construct an orbital at a particular location r, are defined 
as 



F L ^( n ) = 



k 



(3.1) 



where L(S) represents for large (small) component, k de- 
notes number of GTOs, Af correspond to normalization 
factor for each GTO and % is an arbitrary parameter 
which has to be chosen suitably for orbitals from differ- 
ent £ symmetries. To get more flexibility in optimization 
of our basis sets, we use the even tempering condition by 
defining two more parameters £ and v as 



The radial grid points rt are defined as 

n = r [e h ^ - 1], 



(3.2) 



(3.3) 



with ro is the starting radial function taken inside the 
nucleus to be 2 x 10~ 6 at which the wave functions be- 
come finite and h is a step size which is defined by taking 
maximum radial function r max as 150.0 au and total grid 
points 1000. 

We have considered 40 GTOs for each I symmetry or- 
bitals and the considered £ and v are given in Table Q] for 
different I values. Due to limitation over computational 
resources and negligible contributions from the high ly- 
ing virtual orbitals, we have taken up to 24 orbitals from 
s, p, d symmetries and 17 orbitals from /, g symmetries 
to construct active space for RCC calculations. 

Also, the orbitals are generated by accounting the fi- 
nite size of the nucleus assuming a two-parameter Fermi- 
nuclear-charge distribution given by 



Po 



1 + e( ri_c )/ a ' 



(3.4) 



where po is the density for the point nuclei, c and a are 
the half-charge radius and skin thickness of the nucleus. 
These parameters are chosen as 



and 



a = 2.3/4(Zn3) 



^' rms ^ 



(3.5) 



(3.6) 
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where r rms is the root mean square radius of the corre- 
sponding nuclei which is taken as 3.61 fm (27j . 



B. Calculation of atomic wave functions 

To calculate matrix elements of the hyperfine inter- 
action Hamiltonian, we use the RCC method where we 
define atomic wave functions for the considered states 
with valence orbital denoted by v as [28, 29] 



\* v ) = e T {I + 



(3.7) 




(i) 



(iii) 



Hn „ 



V 




iv 



V 



where the DF wave function is constructed as = 
a£|<I>o) with |$ ) is the DF wave function for the closed- 
shell configuration [3p 6 ] in the considered K atom. In 
the above expression, T and S v are the excitation opera- 
tors that accounts for core and core-valence correlations 
to all orders, respectively. Since K is a small size atom, 
hence correlation effects among electrons are expected to 
be less. Therefore role of the higher order configurations 
in determining atomic wave-function could be negligible. 
On the other-hand consideration of these configurations 
are computationally very expensive. Owing to this fact 
we account only all possible single and double configura- 
tion excitations to all orders (known as CCSD method) 
by expressing the above operators in the Fock space rep- 
resentation as 



T = T 1 +T 2 = 



i»i 

b 



a.p 



ab,pq 



(3.8) 



S v = Sl v + S 2v = ^ a l a vS P v + \ X! a i a l a b a vSy, 



a.p 



ab,pq 



(3.9) 



where the (a,b,c - ■ ■), (p,q,r---) and (v) subscripts of 
the second quantized operators represents core (hole), 
particle (virtual) and valance orbitals, respectively. How- 
ever expanding Eq. (|3.7p using these CCSD operators to 
all non-linear terms give rises contributions from higher 
excitations. We determine the above t and s v coefficients 
which correspond to the excitation amplitudes using the 
following equations 



(* L \{H N e T }\$ )=0 (3.10) 



and 



(*5\{H N e T }S v \* v ) = -(*t\{H N e T }\* v ) 

+ (^\S V \^ V )AE V , (3.11) 

with the superscript L(=l,2) representing the single and 
double excited configurations from the corresponding DF 
states, the wide-hat symbol denotes the linked terms, 
AE V is the attachment energy of the valence electron 



FIG. 1: Typical Goldstone diagrams representing leading- 
order triple excitations over the CCSD method. Double 
arrows in the diagrams represents valence electron v, and 
the lines with upward (downward) arrows represents parti- 
cle (hole) orbitals. 



v and Hn denotes the normal ordering atomic Dirac- 
Coulomb Hamiltonian H which is taken as 



H = 5>"rP* + (ft - l)c 2 ] + J2 — > 



(3.12) 



where a and ft are usual Dirac matrices, c is the velocity 
of light. AE V is evaluated by 



AE V = ($ v \{H N e T }{l + S v }\$ v ) 



(3.13) 



To improve the quality of energy and calculation of wave 
functions due to the dominant triple excitations contain- 
ing the valence orbital, we define a perturbation operator 
S3V by contracting Hn with T 2 and S 2v operators as 



s 3v{s vbc ) - 



HnT 2 + HnS 2v 



f-r — Cfc — £c 



(3.14) 



with sffl correspond to excitation amplitudes and is 
the DF energy of the electron in the i th orbital. This 
operator is considered as a part of S v operator in Eq. 
(|3 . 13|) to get additional contribution to AE V . Since AE V 
is involved in Eq. (|3.11|1 . we solve both the equations 
simultaneously in an iterative procedure. This approach 
is generally referred as CCSD(T) method 30]. The dia- 
grammatic representation of these excitations are shown 
in Fig. □ 

The expectation values due to the hyperfine interaction 
operators have been evaluated using our RCC method by 



{T ( k)) = 



(^ V \T^\^ V ) 

(*o|*„) 

{l + Sl}N {l + S v } 

($ v \{i + s\ v + sl}¥* ] {1 + s lv + s 2v }\$ v ) 



{1 + Si + Sl v }N {l + S lv + S 2v } 



(3.15) 
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where T? 5 = (e Tt T e (fe) e T ) and JV = e Tt e T . Gener- 

ally, both r; ' and iV in our RCC approach are non- 
terminating series. These terms are terminated keeping 
terms minimum up to fourth order in perturbation. De- 
scription of this procedure has been given in the previous 
works [3H - I33I ] . Contributions from normalization of the 
wave functions (Norm) are estimated explicitly in the 
following way 

AT rm=<* K |TW|*, ) )| r i_-lJ (3.16) 
where N v = {1 + S\ v + St v }N Q {l + S lv + S 2v }- 

IV. RESULTS AND DISCUSSIONS 

Our aim is to obtain B/Q values more accurately in 
different states of K atom so that they can be combined 
with the available precise experimental results for B to 
estimate Q. In order to verify the accuracies of B/Q 
results from our calculations, it would be felicitous to 
test the accuracies of the wave functions in the nuclear 
region. Owing to the fact that our calculation proce- 
dure deals with many numerical computations at differ- 
ent stages, along with that it correlates with higher exci- 
tations configurations indirectly, hence it would be very 
difficult to estimate uncertainties from the used numer- 
ical methods and approximations taken at the level of 
excitations. With the intention of verifying accuracies of 
the wave functions in the nuclear region, we have calcu- 
lated A for many states in K. Assuming that the anoma- 
lous effects due to different nuclear sizes in all the consid- 
ered isotopes are very small, we evaluate A/gi in 39 K and 
determine A values for the corresponding isotopes using 
their respective gi values. We have used experimental 
values of gi{ 39 K) = 0.2609772, gi( A0 K) = -0.324525 
and gi( 41 K) — 0.1432467 @ to estimate these quantities 
ignoring their uncertainties as they will not meddle the 
results within the reported uncertainties. Both the cal- 
culated and experimental results are compared in Table 
HTl We estimate uncertainties in our calculations by con- 
sidering incompleteness of basis functions, contributions 
from the inactive orbitals in the RCC method, higher 
order excitation levels and from the neglected terms in 
the non-truncative series of Eq. (13. 15)) . The upper limit 
to these uncertainties are given in the parentheses of the 
above table. Clearly, our estimated uncertainties are fair 
enough to compare with the available experimental re- 
sults and our assumption for neglecting the anomalous 
effects for different isotopes seem to be reasonable. 

We also compare our A results for 39 K obtained using 
DF, CCSD and CCSD(T) methods with another recent 
calculations [42| in Table [TTT1 In Ref. [42} Safronova and 
Safronova have also used linearized RCC method with 
singles and doubles approximation (SD method) and in- 
cluding important triples effects for some of the states 
(SDpT method). They find large differences between 



their SD and SDpT results in contrast to our finding 
of small differences between our CCSD and CCSD(T) 
results. However, both the calculations reveal that the 
signs of the A values of the 3D 5 / 2 and 5D 5 m states are 
negative, which were not resolved correctly in the mea- 
surements. Moreover, these calculations indicate that 
correlation effects in the considered atom are substantial 
for which an all order perturbative method like ours is 
suitable to determine wave functions accurately. To our 
knowledge, A values are not known experimentally for 
some of the states in 39 K and close agreement between 
the results from both the calculations in these states will 
be very useful to conduct new measurements in the right 
direction. Also, we have given A values for few excited 
states where neither measurements nor theoretical calcu- 
lations are available. 

Following A results, we present now our calculated 
B/Q results using the DF, CCSD and CCSD(T) methods 
in Table IIVI for the states where precise experimental B 
results for 39 K and 41 K are available. We also estimate 
the uncertainties associated with the results obtained us- 
ing CCSD(T) method and present them in the parenthe- 
ses of the same table. The uncertainties are estimated us- 
ing the procedure as we followed for A. To justify that, 
trends of correlation effects for both the properties be- 
have in a similar manner, we present contributions from 
various RCC terms to both A and B/Q results in Table 
|V]for few important states. It is found that, states where 
angular momentum is larger than half, the contributions 
from the correlation effects are always of similar scale in 
magnitude and the contributions coming from OS2v + cc 
are larger than OS± v + cc. This implies that the core- 
polarization effects are large to estimate B/Q results in 
the considered atom which are accounted for, up to all 
orders through OS2V + cc RCC terms in our calculations. 
There are no other calculated results for B/Q available 
to our knowledge in any of the considered isotopes of K 
to compare with our results. 

Both the calculated results for A and B/Q seem to be 
very accurate, moreover A values are in good agreement 
with the available experimental results. There are also 
several experimental results available for B in 39 K as well 
as 41 K which are given in Table [TV] [H H, El]. 
The most precise values are quoted in bold fonts for the 
respective states in the same table and we combine these 
results with our calculated B/Q values to estimate Qs in 
these two isotopes. We obtain three different values of Q 
in 39 K and two values in 41 K. All these estimated values 
agree with each other in their respective uncertainties, 
but the most precise results which are obtained from the 
4P 3/2 state are 0.0625(17) b and 0.0744(10) b for 39 K 
and 41 K, respectively. Here we have used the following 
expression to evaluate the net uncertainties of Q values 

where we assume C is the extracted value from A/B and 
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TABLE II: Comparison of calculated and available experimental A results in 39-41 K (in MHz). Theoretical As in different 
isotopes are evaluated using the calculated CCSD(T) results of A/gi in 39 K and their respective experimental gi values. 
Uncertainties estimated from the calculations are given in parentheses of our results. 



State 




This work 






Experiment 






39 K 


40 K 


4i K 


39 K 




41 K 


4Si/a 


229.6(2.0) 


-285.5(2.5) 


126.0(1.1) 


230.8598601(3) [3] 


-285.7308(24) [3] 


127.0069352(6) [3] 


iPl/2 


27.4(5) 


-34.04(62) 


15.02(27) 


27.775(42) [2] 
28.85(30) [3] 

27.80(15) [34] 
27.5(4) [35] 

28.859(15) [36] 


-34.523(25) [2] 
-34.49(11) [34] 


15.245(42) [2] 
15.19(21) [34] 
15.1(8) [35] 


4P 3 /2 


5.9(2) 


-7.35(37) 


3.25(16) 


6.093(25) [2J 
6.06(8) [3] 
6.00(10) [38] 
6.13(5) [3_7J 


-7.585(10) [2J 
-7.48(6) [34] 
-7.59(6) [39] 


3.363(25)[2J 
3.40(8) [37] 
3.325(15) [40] 


3£»3/2 


1.0(2) 


-1.25(25) 


0.55(11) 


0.96(4) [1] 


1.07(2) [1] 


0.55(3) [1] 


3£>5/2 


-0.57(5) 


0.711(62) 


-0.314(27) 


0.62(4) [1] 


0.71(4)[1] 


0.40(2) [1] 


4£>3/2 


0.686(4) 


-0.853(5) 


0.377(2) 








4^5 /2 


-0.332(2) 


0.413(2) 


-0.182(1) 








5Si/2 


55.0(1.0) 


-68.4(1.2) 


30.18(55) 


55.50(60)[3J 






5Pl/2 


8.8(5) 


-10.98(62) 


4.84(27) 


9.02(17) [3] 






5P3/2 


1.9(2) 


-2.37(25) 


1.04(11) 


1.969(13) [3] 
1.95(5) [38] 


-2.45(2) [3] 


1.08(2) [3] 


5£>3/2 


0.39(5) 


-0.489(62) 


0.22(27) 


0.44(10) [3] 






5D 5/2 


-0.171(7) 


0.213(9) 


-0.094(4) 


±0.24(7)[3J 






6S1/2 


21.1(6) 


-26.24(75) 


11.58(33) 


21.81(18)[3] 




12.03(40) [3] 


6P1/2 


3.9(4) 


-4.87(50) 


2.15(22) 


4.05(7) [3] 






6P3/2 


0.9(2) 


-1.17(25) 


0.51(11) 


0.886(8) [3J 






6-D3/2 


0.24(5) 


-0.297(62) 


0.131(27) 


0.25(3) [41] 
± 0.2(2) [3] 






6^5/2 


-0.112(5) 


0.139(6) 


-0.061(3) 


-0.12(4) [41] 
± 0.10(10) [3] 






7Si/ 2 


10.3(5) 


-12.83(62) 


5.66(27) 


10.79(5) [3] 






7Pl/2 


2.1(3) 


-2.61(37) 


1.15(16) 


± 2.18(5)[3J 






7P3/2 


0.50(3) 


-0.62(37) 


0.28(16) 








7D 3/2 


0.15(1) 


-0.19(12) 


0.083(5) 








7D 5/2 


-0.068(5) 


0.085(6) 


-0.037(3) 








8S1/2 


5.8(2) 


-7.27(25) 


3.21(11) 


5.99(8)[3J 






8P1/2 


1.3(1) 


-1.56(12) 


0.69(55) 








e p 

»-f3/2 


U.oUl^o J 


n Q7/i ( «\ 
— U.O ( 4(^0 J 


U.lDO^o ) 








8-D3/2 


0.101(2) 


-0.126(2) 


0.055(11) 








8A/2 


-0.040(1) 


0.050(1) 


-0.022(1) 








9Sl/2 


3.4(3) 


-4.25(37) 


1.88(16) 










1.0(1) 


-1.19(12) 


0.53(55) 








9P 3 /2 


0.230(4) 


-0.286(5) 


0.126(2) 








9^3/2 


0.334(4) 


-0.415(5) 


0.183(2) 








9£>5/2 


-0.148(2) 


0.184(2) 


-0.081(1) 








105 1/2 


2.2(3) 


-2.68(37) 


1.18(16) 


2.41(5) [3] 







and 

= 1.2173 ±0.0001. (4.3) 

Using the measured and Q value of 41 K, we get 

a new Q value for 39 K as 0.0612(8) b. Considering both 
the values of Q in 39 K, we restrict the lower and upper 
limits of Q( 39 K) to 0.0608 b and 0.0620 b, respectively. 



5A, 5B and SC are their respective uncertainties. 

Among both the new Q values in 39 K and 41 K, the 
relative uncertainty in Q of 41 K is small. Moreover there 
are also experimental results for the ratios of Q values 
between 39 K, 40 K and 41 K are available as 144, 46, 47] 



= -1-244 ±0.002 (4.2) 
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TABLE III: Comparison between different theoretical results of A in 39 K (in MHz). 



State 




Others [42] 






This Work 




DF 


SD 


SDpT 


DF 


CCSD 


CCSD(T) 


/ 


146.91 


237.40 


228.57 


146.794 


229.573 


229.556 


4Pl/2 


16.616 


28.689 


27.662 


16.616 


27.247 


27.371 


4P 3 /2 


3.233 


6.213 


5.989 


3.234 


5.886 


5.913 


30 3 /2 


0.447 


0.983 


1.111 


0.447 


1.006 


1.003 


3£>5/2 


0.192 


-0.535 


-0.639 


0.192 


-0.574 


-0.572 


4£>3/2 


0.281 


0.678 




0.281 


0.690 


0.686 


40>s/2 


0.120 


-0.307 




0.120 


-0.334 


-0.332 


5<Si/2 


38.877 


56.102 


54.817 


38.847 


55.070 


54.981 


5Pl/2 


5.735 


9.202 


8.949 


5.735 


8.755 


8.827 


5P 3 /2 


1.117 


1.988 


1.932 


1.117 


1.887 


1.903 


503/2 


0.168 


0.409 




0.168 


0.396 


0.393 


5^5/2 


0.072 


-0.167 




0.072 


-0.173 


-0.171 


6S1/2 


15.759 


22.025 


21.609 


15.105 


21.167 


21.10 


6P1/2 


2.629 


4.066 


4.014 


2.629 


3.874 


3.918 


6P 3 /2 


0.512 


0.874 


0.866 


0.512 


0.928 


0.937 


60>3/2 


0.105 


0.253 




0.104 


0.241 


0.239 


60>5/2 


0.0448 


-0.0975 




0.045 


-0.113 


-0.112 


7iJl/2 


7.900 


10.876 


10.690 


7.894 


10.363 


10.317 


7Pl/ 2 


1.417 


2.191 


2.140 


1.417 


2.066 


2.095 


7^3/2 


0.276 


0.473 


0.462 


0.276 


0.495 


0.502 


7D 3/2 


0.0685 


0.1644 




0.067 


0.154 


0.152 


7^5/2 


0.0293 


-0.0611 




0.028 


-0.069 


-0.068 


8<S'i/2 


4.511 


6.156 


6.057 


4.536 


5.880 


5.847 


8-P1/2 








0.855 


1.236 


1.257 


8-P 3 /2 








0.166 


0.296 


0.301 


80*3/2 








0.048 


0.102 


0.101 


SD 5/2 








0.019 


-0.040 


-0.040 


95 1/2 


2.814 


3.818 


3.759 


2.685 


3.444 


3.420 


9Pl/2 








0.695 


0.945 


0.958 


9P 3 /2 








0.136 


0.227 


0.230 


9D 3/2 








0.189 


0.338 


0.334 


9A/2 








0.083 


-0.151 


-0.148 


105 1/2 


1.871 


2.529 


2.491 


1.878 


2.171 


2.154 



TABLE IV: Calculated B/Q values (in MHz/fc) from our DF, CCSD and CCSD(T) methods and accurately known experimental 
B results (in MHz) in 39 K and 41 K. The most precise B results are given in bold fonts and Q values are estimated by combining 
these results with the calculated B/Q values in both the isotopes. 



State 


DF 


CCSD 


Theory 

CCSD(T) 


Experiment 
B( 39 K) 


qC 9 k) 


Experiment 
B( 41 K) 


Q( 41 K) 


4P 3 /2 


23.000 


44.392 


44.6(5) 


2.786(71) [2] 


0.0625(17) 


3.351(71) [2] 












2.72(12) [36] 




3.34(24) [37] 












2.83(13) [43J 




3.320(23) [40J 


0.0744(10) 










2.9(2) [38J 








5P 3 /2 


7.866 


13.833 


13.9(4) 


0.870(22) [3] 


0.0626(22) 


1.06(4) [37] 


0.0763(36) 










0.92(10) [38] 








6P3/2 


3.340 


6.285 


6.4(5) 


0.370(15) [3] 


0.0578(51) 







Therefore, we recommend Q value of 39 K as 0.0614(6) b. 
Now with this most precise Q value and the above ratios 
of Q values between different isotopes, we get precise Q 
values for 40 K and 41 K as -0.0764(8) b and 0.0747(7) 
b, respectively. There are also other reported Q values 
which we have compared them with ours in Table I VII 
As seen from the table, there are three other works re- 



port these results [ll|, |lj, |44| . Apart from Ref. [4||, the 
calculations carried out in these works are rigorous and 
results reported in are the latest. Q values reported 
in 11( matches with our estimated values with some over- 



laps within the predicted uncertainties, however results 
reported in [12| disagree with us. In both these theoreti- 
cal works, they have determined electric field gradients at 
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TABLE V: Contributions from different RCC terms to A (in MHz) and B/Q (in MHz/6), whereas c.c. stands for complex 
conjugation. 



State 


O 




OSiu - 


h c.c. 


os 2 „ 


+ c.c. 


Others 


Norm 








B/Q 




B/Q 


A 


B/Q 


A 


B/Q 


A 


B/Q 


451/2 


146.576 




46.148 




33.0457 




6.707 




-2.920 






16.341 




5.500 




4.820 




n aaa 
0.9UU 




-0.190 




4P 3 /2 


3.186 


22.662 


1.062 


7.562 


1.430 


13.386 


U.z/O 


1 OO A 

1.280 


-0.040 


- 0.305 


3I>3/2 


0.464 


1.168 


0.430 


1.080 


-0.059 


2.947 


n lor 


A AOA 


-0.017 


- 0.090 


3£» 5 /2 


0.199 


1.666 


0.183 


1.527 


-0.766 


4.223 


— U.l la 


A 1 0^7 
0.12/ 


0.010 


- 0.128 


5<Si/2 


38.770 




7.898 




8.200 




0.725 




-0.612 




5*1/2 


5.683 




1.535 




1.571 




0.100 




-0.062 




0-T3/2 


1 nsQ 


7 7^9 


fl 9Qfi 
u.zyu 


Z . 1 11 




o.yu 1 


0.068 


0.206 


n m ^ 

— U.UIO 


n nofi 

u.uyu 


6Sl/2 


15.043 




2.882 




3.150 




0.246 




-0.217 




6*1/2 


2.615 




0.633 




0.630 




0.067 




-0.027 




6*3/2 


0.510 


3.631 


0.122 


0.872 


0.276 


1.825 


0.037 


0.070 


-0.006 


-0.043 


75i/ 2 


7.902 




1.148 




1.309 




0.062 




-0.104 




8S1/2 


4.540 




0.592 




0.746 




0.026 




-0.057 




9Si/ 2 


2.687 




0.316 




0.438 




0.011 




-0.032 




10Si /a 


1.880 




0.013 




0.291 




-0.014 




-0.015 





TABLE VI: Reported values of nuclear quadrupole moments 
Q in b for 39 K, 40 K and 41 K from various studies. 



Isotope 



This work 



Others 



K 



°K 



K 



0.0614(6) 

-0.0764(8) 
0.0747(7) 



0.0585 [12] 
0.0601(15) [11] 

0.049(4) [44] 

-0.073 a [12] 
-0.0749(19) [11] 

0.0711 [12] 
0.0733(18) [11] 



"Accuracy is expected to be better than 1 percent 



the nucleus to extract the nuclear quadrupole moments 
and the results are model independent. But both the 
calculations are less rigorous than the present calcula- 
tions. In Ref. |ll| . Sundholm and Olsen have used a non- 
relativistic large scale finite-element multi-configuration 
Hartree-Fock configuration interaction (MCHF) method. 
The core contributions were estimated from the core- 
valence correlation calculations and the relativistic cor- 
rections are accounted separately from the DF calcula- 
tion. Contrast to this work, we have considered the core 
correlation and core- valence correlations to all orders and 
relativistic effects are included to all orders through the 
RCC method. In fact, their truncative CI method is 
known to have size-consistent problem |48j against our 
RCC method. On the other-hand, level of approxima- 
tions employed to carry out calculations in Ref. [12[ by 
Kello and Sadlej are comparable to the present work. In 
their work, calculations are performed with a scalar rel- 
ativistic Hamiltonian using Douglas-Kroll approach and 
CCSD(T) method is used to account the correlation ef- 



fects. Since this approach is better than the above MHCF 
method, the results reported in 12[ were considered to be 
more accurate and estimated to be within 1% accuracy. 
However, we find the reported values in [l2j are smaller 
in magnitudes as compared to our estimations. 

It will be interesting to see further theoretical stud- 
ies of B/Q results to draw comparison with the present 
work. Further, we combine our calculated B/Q results 
with the newly obtained Q values to determine theoreti- 
cal results for B in many states of the considered isotopes 
of K. These results are given in Table lylll We have also 
neglected here the anomalous effects in calculated B/Q 
values for different isotopes due to their negligible roles. 
If some of these B results can be measured more precisely 
than the reported results, then combining those results 
with our calculated B/Q values will definitely give rise 
better Q values in these isotopes. In fact, more precise 
theoretical calculations of B/Q results in these isotopes 
can also give rise to more accurate Q values in K. 

As can be noticed in Table Ivnl our estimated B results 
for different states in K isotopes agree very well in most 
of the states except with the 3D states. We suggest to 
carry out further measurements of B in these states to 
ascertain our results. Theoretical results for B are also 
given in many excited states for the first time which can 
be verified by the future measurements. 



V. CONCLUSION 

We have employed relativistic coupled-cluster method 
to calculate matrix elements of the hyperfine interaction 
Hamiltonians in potassium atom. By performing calcula- 
tions of the magnetic dipole hyperfine structure constants 
in this atom, we have tested the accuracies of the wave 



8 



TABLE VII: Comparison of estimated and experimental B results in 39 41 K (in MHz). 



State 


(s/Q) (Tft) 




This Work 






Experiments 








39 K 






39 K 


40 K 




4P 3 /2 


44.6(5) 


2.738(41) 


-3.41(5) 


3.332(49) 


2.786(71)[2J 
2.9(2) [38] 
2.72(12) [37] 
2.83(13) [3] 


-3.445(90) [2] 
-3.23(50) [34] 
-3.5(5) [39] 


3.351(71)[2J 
3.34(24) [3JJ 
3.320(23) [40J 


3£>3/2 


5.2(5) 


0.319(31) 


-0.40(4) 


0.388(38) 


0.37(8) [1] 


0.4(1) [1] 


0.51(8)11] 


3£>5/2 


7.4(4) 


0.454(25) 


-0.565(31) 


0.553(30) 


<0.3 [1] 


0.8(8)[1] 


< 0.2 [1] 


4£>3/2 


2.35(4) 


0.144(3) 


-0.180(4) 


0.176(3) 








4£> 5/2 


3.35(4) 


0.206(3) 


-0.256(4) 


0.250(4) 








5P 3 /2 


13.9(4) 


0.853(26) 


-1.06(3) 


1.038(31) 


0.870(22) [37, 38] 
0.92(10)[38J 


-1.16(22)[39J 


1.06(4) [37] 


5£>3/2 


1.16(8) 


0.071(5) 


-0.088(6) 


0.087(6) 








5£> 5 /2 


1.65(10) 


0.101(6) 


-0.126(8) 


0.123(8) 








6P3/2 


6.4(5) 


0.393(31) 


-0.489(39) 


0.478(38) 


0.370(15) [45] 






6-D3/2 


0.72(5) 


0.044(3) 


-0.055(4) 


0.054(4) 


0.05(2) [41] 






6D5/2 


1.02(6) 


0.063(4) 


-0.078(5) 


0.076(5) 








7P 3 /2 


3.4(3) 


0.209(19) 


-0.260(23) 


0.254(23) 








7-D3/2 


0.44(3) 


0.027(2) 


-0.034(2) 


0.033(2) 








7D 5/2 


0.62(4) 


0.038(2) 


—0.047(3) 


A C\ A a { \ 

0.046(3) 








8P3/2 


2.0(2) 


0.123(12) 


-0.153(15) 


0.149(15) 








8^3/2 


0.29(2) 


0.018(1) 


-0.022(2) 


0.022(2) 








8C5/2 


0.36(2) 


0.022(1) 


-0.028(2) 


0.027(2) 








9P 3 /2 


1.5(2) 


0.092(12) 


-0.115(15) 


0.112(15) 








9£>3/2 


0.92(6) 


0.056(4) 


-0.070(5) 


0.069(5) 








9D 5/2 


1.35(8) 


0.083(5) 


-0.103(6) 


0.101(6) 









functions in the nuclear region. These wave functions 
were further used for the electric quadrupole hyperfinc 
interaction studies. By combining our calculations with 
the corresponding measurements, we obtained the nu- 
clear quadrupole moments as 0.0614(6) b, —0.0764(8) b 
and 0.0747(7) b for 39 K, 40 K and 41 K, respectively. These 
results agree with one of the previous work but do not 
agree with others including the latest reported results. 
After obtaining nuclear quadrupole moments, we substi- 
tuted them to obtain electric quadrupole hyperfine struc- 
ture constants in many states and found very good agree- 



ment with the experimental results except for the 3-D 
states. Also, we have given some of the results that were 
not reported earlier. We suggest further studies of the 
considered properties to ascertain our results. 
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